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FOREWORD 


The Space and Informations Systems Division (SID) of North American 
Aviation, Inc. (NAA) under Contract AF 30(602)-3628 with the Rome Air 
Development Center (RADC) of the United States Air Force agreed to perform 
@ ten month study designed to develop digital computer techniques in two 
areas of iaterest to the RADC tracking facility. First, a differential 
correction geocentric orbit computation program was to be prepared for 
reducing observed data which will operate in a near optimum manner at 
the RADC computer center. And second, e@ computational logic which could 
be utilized in the tracking process for driving the tracking antennae in 
an open loop mode was to be prepared. This second program would employ 
general perturbations theory in the definition of the predicted trajectory. 


This report was prepared as partial documentation of the first task. 

The contents present the prog.cam logic and FORTRAN listings for adata handling 
preliminary processor to be utilized in conjunction with the differential 
corrections program (DCP) described in a previous volume (SID 65-1203-1). 
This preprocessor is designed to accomplish two functions, one of which 
is required in the main program. The first function involves the sorting 
of the data and the arrangement in a chronological format. This operation 
provides a convenient means of incorporating data from several stations 
and assures that recursive form of the data filter employed in the DCP 
will perform as desired. The second functim provided is a preliminary 
smoothing and editing of the raw data. This operation is accomplished by 
fitting the raw data acquired over a short interval of time to a parabolic 

xe in the sense of least squares for the purpose of eliminating random 
scatter. The smoothed data are edited by recording only the mid point of 
the segment on a magnetic tape along with all of the information necessary 
to identify it. This latter step is performed in the interest of 
computational efficiency with particular thought given to optimizing the 
program for the RADC facility. 


This contract has been w:anaged at NAA S&ID by Mr. J. A. Hill and 
directed by Mr. G. BE. Townsend. Mr. Townsend also conceptually designed 
the rationale for the preprocessor and assisted in its preparation. The 
formal logic for the program in addition to its development and checkout 
are the products of Mr. C. C. DeBilzan. 


The assistance offered by RADC personnel under the direction of 
Mr. Gordon Negus (Program Manager) is gratefully acknowledged. 
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ABSTRACT 


The Preprocessor for the Differential Corrections Geocentric Orbit 
Computations Program (FS4~305A) is a data handling routine required in the 
reduction of satellite observations, and the computation of the associated 
trajectories, This program reads a specially generated magnetic tape con- 
taining coordinate data observations fiom a series of tracking stations. The 
data is smoothed, using least squares criteria, and then ordered. Primary 
ordering is time and secondary ordering is by station. Output includes a 
magnetic tape which is required as input for the Differential Corrections 
Geocentric Orbit Computations Program (FS4-305). 
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INTRODUCTION 


FS4-305A is a special purpose, FURTRAN IV, IBM 7094 program which 
was written and checked using the standard North American Aviation monitor 
system (NAASYS-version 13). The primary function of this program is to 
act as the interface between the tracking stations and the Differential 
Corrections Geocentric Orbit Computations Program (FS4-305). As illustrated 
in figure 1, this function is executed by performing a preliminary smoothing 
of the data to eliminate random scatter and reduce the magnitude of the 
processing task (by eliminating some of the data) and by arranging the 
data in the format expected by FS4-305. The arrangement involves the 
construction of a seven word record for each smoothed set of observations, 
identifying the time, the station, the type of data and the observation 
vector. This data is presented on a magnetic tape, chronologically. 


The data will be read from this tape in the differential 
corrections prograi (DCP) described in SID 65-1201-1 by subroutine DATAPE. 
This routine is the sole link between the preprocessor and the DCP, and as 
such it assures complete consistency between the two programs. Accordingly, 
its design and operation are discussed in this document and in the documen- 
tation of the DCP. 


Tt was anticipated that the amount of data being processed may, at 
times, require storage capabilities exceeding core. Further, it was 
recognized that the dimensions of certain key arrays within the program 
are dependent upon the type of data being processed. Therefore, provision 
has been made for auxiliary tape storage (mechanized by the program only 
if required) and variable dimensioning has been utilized. Proper adjustment 
of array dimensions, either through knowledge of certain key parameters 
beforehand or with the aid of diagnostics written under program control, 
assures that the auxiliary tape storage mode will usually not generally 
require mechanization. 


Though the program was developed on the IBM 7094 itwas also designed to be 
operated on CDC equipment (specifically the CDC 1604) and other systems 
possessing FORTRAN capabilities but not utilizing NAASYS. Consequently, to 
minimize possible system incompatabilities, a “straight forward” approach 
wap utilized in both setting up the program logic and coding. Although this 
philosophy preciuded the use of machine language routines, it is beileved that 
the programs efficiency has not been appreciably impaired. Further, operation 
on new systems such as the IBM 360 is assured. 
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WHIGH SERVES AS INTER-FACE 


FIGURE ! 


STEPS REQUIRED TO TRANSMIT DATA FROM TRACKING 
STATIONS TO FS 4-305, {OIFFEREN TIAL CORRECTIONS GE@CENTRIC 
ORBIT COMPUTATIONS PROGRAM), 


®READS INPUT DATA FROM 
VARIOUS STATIONS, 
‘OUTPUTS ALL BATA ON 
BETWEEN TRACKING STATIONS]| A SINGLE MAGNETIC TAPE 
AND F6 4-305 A (FORTRAN IS, BINARY MODE) 
WITH A S°ECIFIED FORMAT. 


SINA 
RAW DATA 


SPECIFIC PURPOSE ROUTING 


eee 


FILE PER STATION PER PASS 


RAW DATA 
(FORTRAN IW, BINARY MODs) 


(MAG) 


RAW DATA ARRANGED ONE { 


* READS RAW DATA FILE 
BY FILE. 
* SMOGTHS RAW DATA FILE 
BY FILE, 
* COMBINES SMOOTH DATA 
FILES ANO ORDERS DATA, 
* MAKES NECESSARY UNIT 
CONVERSIONS. 
* WRITES THE SMOCTHED AND 
OBERED SATA ON TAPE TO 
GE READ BY MAIN PROGRAM, 


*STN ORDER 
*CONVERSION 


FS4-305A COEF. 


PRE-PROCESSER FOR 
DIFFERENTIAL CORRECTIONS 
PROGRAM. 


T RAW AND/OR 
SMOOTHED 
DATA. 


SMOOTHED 
ANC 
ORDERED 
DATA 
(MAG) 


SMOOTHED AND ORDERED 
DATA ARRANGED WITH A 
SPECIFIES NUMBER OF 
POINTS PER LOGICAL RECORD 
(FORTRAN TZ, BINARY MODE.) 


* SUBROUTINE OATAPE READS 
SMOOTHED AND ORDERED 
DATA POINT BY POINT. 


FSY~- 305 
DIFFERENTIAL CORRECTIONS 
GEOCENTRIC ORBIT 

COMPUTATIONS PROGRAM. 
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SUMMARY OF PROGRAM CHARACTERISTICS 


This data handling routine has been casigned to assure complete 
consistancy between the differential corrections program and the sources 
of the observation data. For this reason, before proceding 
with the discussion of the logical structure of the routine, it is singularly 
instructive to outline the requirements imposed on the output of the 
preprocessor and some of the features which are desired of it during its 
operation. An itemization of such features is presented below: 


@Output is a magnetic tape required as input to FS4-305. 


e Output tape contains a series of points, each point consisting of seven 
words indicating time, station, data type, and up to three coordinates, 


e Input requirements include a specially generated magnetic tape containing 
coordinate data from tracking stations recorded at about one second intervals.* 
a, coordinate data may be any combination of the following types. 
1. Range 
2. Range rate 
3. Azimuth, elevation 
4k, Range, range rate 
5. Range, azimuth, elevation 
6. Range rate, azimuth, elevation 
eUp to ten separate tracking stations permitted (single satellite). 


eOther than exceeding data storage capabilities, no limiu on number of passes 
per station. 


eSmooths data using least squares criteria. 
eOrders data. Primary ordering is time, secondary ordering is station. 


eObservational sphere of influence of ad‘acent stations may overlap. 
i.e. More than one station may record observations at the same time. 


eFORT! -N IV (IBSYS-version 5, NAASYS-version 13). 
eVariable dimensions. 

eCOMMOGN not used. 

eRaw and/or smooth data print options. 

eAuxiliery tape storage mechanized by program if required. 


* If 4? >> / subroutine FIT should be reviewed and fewer than 20 data points 
smcotned at a time, ~3- 
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@Requires special tapes be mounted. 


eRequired logical magnetic tape units: 
tape 5, System input 
tape 6, System output 
tgpe 8, System open. Used for inputting raw data observations from 
vy tracking stations. 
tape 9, System open. Used for: 
1.) temporary data storage, 
2.) inputting smoothed and ordered data. 

o Due to the general nature and flexability of the program, the maximum 
number of raw data points that may be processed is difficult to estimate. 
Generally, the upper limit will be reached if the number of smoothed 
data points exceeds the dimension of the STN array. 


e Requires approximately 13465, 3% 594119 locations excluding locations 
required for STN,A,AA arrays end I/¢ buffers. 


« Loading and execution time for the sample problem was less than 2 minutes 
on the IBM 7094, 
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METHOD 


FS4-305A reads a specially generated magnetic tape containing 


observations of a Satellite from a series of tracking stations. The . 
deta is smoothed using least squares criteria and chronologically ordered. 


After making the required unit transformations, the data is output on 
magnetic tape in the format requin2d by the main program, FS4-305. 
Provision is made for optional printout of the raw and/or smoothed data. 


This function is executed by mechanizing the following routines: 
PRESET presets input data 

REED reads input data 

NBLANK tests for blank input deta fields 


PRINTI prints input data 


PROCES program driver, storage allocation, tape mechanization 
NPT number of rew data points per file 

PRINTR prints raw data 

CHRONT sequences data chronologically 

FIT sets up data for smoothing 


SMOOTH smooths data 
MATMPY matrix multiply 


MAT INV matrix inversion 


CHOPSE check for singular: matrix 
TRANS P matrix transposition 
GSPRT orders data 


CHANGE unit conversion 
TWIF converts time to days 


BESSEL refers time to the epcch of 1950.0 


ANGMgD insures angeles are positive,.mod 360 degrees 
CXA converts degrees to radians 


CXAPLB converts doppler reading to range rate 
PRINTS prints smooth data 


ab 
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The following pages outline noxenclature, method and other pertinent 
details for each subroutine. An apereciation for the complete program logic 
may be gained by referring to the {ow diagrams of the program driver 
routine PROCES. Details concerning the required input tape may be found 
in the section concerning program vyeration whereas the cutput tape is 
described in the program output seciion. A description of the routine 
required te read the output tape (SUBROUTINE DATAPE) is enclosed as appendix 
I. (This routine is the interface between FS4-305A and FS4~305. It is an 
integral portion of F54-305). 
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MAIN routine 


Purpose: Main routine 

Deck Name: MAIN 

Subroutines Required: PRESET (presets input data) 
RSED (reads input data) 
PRINTI (vrints input data) 
PRECES (program driver) 

Functions Required: None 


excluding locations required for STN,a, 


Approximate Deck Length: 166, 3 118 
and AA arrays. 


10 


Method: 


After setting the variable dimensions, input data is initialized, read 
in, and printed. Controi is then transferred to SUBR@UTINE PR@CES which 
serves as the driver program, 


See the sample problem and "Program Operation" section for assistance in 
setting the variable dimensions MAXSTN, MAXA, and MAXAA. 
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Subroutine PRESET 


Purpose! Presets input data and certain program constants, 


Deck Name: PRSET 
Calling Sequence: CALL PRESET (PR, PS, STNAME, C, XJDREF) 
Input/Output: 


FORTRAN | 


Name Description Preset Value 


Raw data print flag 
non zero, print 
zero, no print 


Smooth data print 
Lag 

non zero, print 

zero, no print 


STNAME Station sorting 
order 


Conversion 
coefficients 


XJDREF 


Subroutines Required: 


c(l, J) 
c(2, J)} 
C(3, J) 
ctu, J) 


Pregram reference 
Julian date (zero 
hour U.T. ) (this 

reference will be 


altered in BESSEL to 


refer time to the 


epoch of 1950.0, J.D.| 


2433282,423) 


None 


dates for the sake of numerical significance, 


a1 


Range-no conversion 
ange rate - 

No conversion 

Azimuth, elevation 

degrees to radians 


(33281,.5) * 


* Note that the leading characters "24" have been omitted from all Julian 
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Functions Required: None 


Approximate Deck Length: 122. = 8255 
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Purpose: 


Deck Name: 
Calling Sequence: 
Input /Output: 


FORTRAN 
Name 


STNAME 


Subroutine REED 


Special purpose input routine. As applied to this 
Program, subroutine REED performs the function of 
inputting conversion coafficients and numbers 
required in the ordering of the raw data files (see 
"Program Operation"), These operations are essen- 
tial to insure compatability between this program 
and the Differential Corrections Program, 


RD 


CALL REED (STNAME, C, PR, PS) 


STNAME(J) indicates the sorting order of 
the J-th station. Ordering must be con- 
sistent with the ordering of stations in 
the differential corrections program. 


Coefficients required for unit conversions 
of the observed data. Units must be com- 
patible with those employed in the 
differential corrections program (Km, 
Km/sec, rad) 


Raw data print indicator 


Smooth data print indicator 


Subroutine Required: None 

Functions Required: NBLANK (tests for blank input data) 
Approximate Deck Length: 222, = 148), 

Restrictions: It has been assumed that a blank field is 


assigned the value -0,0, Consequently, 
when inputting data, a zero should be des- 
ignated 0 e 8) 


wipe 
SID 65-1203- 2 


Method: 


Prior to entering, all. input data has been preset at nominal values by 
subroutine PRESET, This routine provides the facility for reading a variable 
number of pieces cf data into the program, replacing the preset values. If 
the preset value is desired, the corresponding input data field is left blank. 
Specific details concerning the input format and instructions may be found in 
the "Program Operation" section, 
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SUBROUTINE REED (STNAME , C,PR, PS) 


READ PRINT 
COA/TROLS 


DO FOR £ACH 
DATA CARO 


READ A OATA CARO 


READ STATION ORDER 
TEST FOR BLANK DATA 
FIELO. IF BLAKIK, OO 
NOT READ DATA /AMTO 
CORE. 


N&: N&LANALDE)? 
$0 

N8:=0? 
SS 


yéS 


ae | Eee! | LJ pd ee 


READ CONVERSION 
COEFFICIENT. 


a a a ae re a Sie | 


SA 


ove | 
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Function NBLANK 


Purpose: Tests for blank input data. 
Deck Name: NBLNK 

Calling Sequence: NB=NBLANK (D) 

Input/Output: 


FORTRAN 
Name Dimensions Description 


Variable to be tested 
blank indicator: 


0, D was a blank 
4, D0 was not a blank 


Subroutine Required: None 
Functions Required: None 
Approximate Deck Length: 50, = 4034 


Restrictions: System must assign the value -0.0 to a blank input 
data field, Zero should be input as 0.0. 
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Subroutine PRINTI 


Purpose: Prints input data 

Deck Name: PRINTI 

Calling Sequence: CALL PRINTI (STNAME, C, PR, PS) 

Input/Output : 

i 
Name Dimensions Description 


STNAME STNAME (J) indicates the sorting order of 
the J-th station. Ordering must be 
consistent with the ordering of stations 


in the differential corrections program. 
Coefficients required for unit conversions 


Raw data print indicator 


Smooth data print indicator 


Subroutines Required: None 
Functions Required: None 
Approximate Deck Length: 306g = 19816 
7 -26- 
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Purpose: 


Deck Name: 


Calling Sequence: 


Inout /Vutput: 


FORTRAN a 


Name 


ag MNAXSTY 


I | MAXA 


NAXAA 


Subroutine PR@CES 


Serves as the program driver routine, 2llocates the 
various data storage requirements and mechanizes the 
required tape drives. 


PRS 


CALL PRACES(STNAMS,MAXSTN , JAXA ,MAXAA,STN,A,AA,XJDR=F, 
PR, PS, C) 


Dimensions Descrintion 


10 StNaK( J) indicites the order in 
which the J-th station is to be 
sorted. 

Dimension of STN array 


Dimension of A array 


Dimension of AA array 
6, MAXSTN STN(I,J) refers to the I-th 
element (coordinate) of the 
J-th point. Used for storage of: 
a.) raw data points (single file 
b.) smoothed points (all files) 
ir auxiliary tape storage 
is mechanized, 
c.) smoothed and sorted points 
(all files) if auxiliary 
tape storage is required, 


6, MAXA A(I,J) refers to the I-th element 
coordinate) of the J-th point. 
Used for storage of: 


\ wad ntea 
Bes SBm00 thed Be ee mvs 


b.) smoothed ind sorted points 
(all files) if auxillary 
storage is not mechanized. 
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FORTRAN 
1/0 Name Dimensions Description 


4, MAXAA AA(I, J) refers to the I-th element 


XJDREF 


Subroutines Required: 


Functions Required: 


Approximate Deck Length: 


Error Indicators: 


Required Logical Magnetic 
Tape Units; 


(coordinate) of the J-th point, Used for 
storage of a single file of smoothed 
points, 


Julian date used for program reference. 
Times associated with the output data are 
in days from XJDREF, See subroutine 
BESSEL for Besselian correction, Not to 
be confused with XDREF, a Julian date 
used as a reference when data is sequenced 
chronologically, 


Raw data print indicator, 
Smooth data print indicator, 


Coefficients required for unit conversions 


PRINTR (prints raw data) 

CHRGNT (sequences data chronologicaily) 
FIT (smooths raw data) 

GSORT (orders smooth data) 

CHANGE (converts units} 

PRINTS (prints smoothed and sorted data) 


NGPT (determines number of points per raw data 
file) 
2454, = 1324, 4 


‘uymerous error tests and diagnostics are incorporated 
within this routine, 


Tape 5, System input 


Tape 6, System output 


Tape 8, System open. Used for inputting raw data 
observations from tracking stations. 


e 9, System n; Used for: 


ape 
temporary data storage (if required) 
outputting smoothed and sorted data 


n 
r 
a 
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Method: 


Prior to entering Subroutine PROCESS, all required input data has been 
read into the program. A specially generated magnetic tape containing coordi- 
nate data from the various tracking stations has been mounted on logical tape 
unit 8, and logical tape wnit 9 is prepared to receive the program output. 
Specific details concerning the format of these tapes may be found in the 
"Program Operation" section. 


A single physical file (i.e., data from one pass of the satellite as 
observed by any one station) of raw data is read from logical tape unit 8 
into the STN array and, after passing through a bookkeeping region, is 
chronologically ordered. Control is then transferred to Subroutine FIT 
where the data is smcothed in an unweighted, least squares mode. The 
smoothed data, corresponding to the raw data file contained within STN, 
is returned from the smoothing routines in the AA array, This array (AA) 
is reused for temporary storage of the smoothed data for each file being 
processed. Finallv, the data within AA is transferred to the A array and 
the procedure is repeated for each file being processed. 


After all raw data files have been read into the program and smoothed, 
the data is then ordered using either of the following methods. If the 
auxiliary tape storage mode were mechanized, the smoothed data has been 
stored on logical tape 9. This data is read into the STN array and ordered 
by Subroutine GSORT. Primary ordering is time, secondary ordering is by 
tracking station. If the auxiliary storage mode were not mechanized, the 
data to be ordered is in the A array. 


After the data has been smoothed and ordered, it is partitioned into 
logical output tape record size, necessary unit conversions are performed, 
and the data is output (binary mode) on tape 9. The O array corresponds 
to a logical output data record. 
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Function N@PT 


Purpose: Determines the number of points per raw data file. 
Deck Name: NPT 

4 Calling Sequence: NP=N@PtT( IFI , NNUM) 

3 Input/Output: 


Dimensions 


Number of coordinates, per point, 
excluding time. 


Total number of words per raw 
data file. 


Number of data points per file. 


0 Sate — rarwrnieieiBerer ane 


Subroutines Required: None 

Functions Required: None 

Approximate Deck Length: 234, = 156) 

Srror Indicators: If NP is non-integer, an error message 
is printed and execution is terminated 
with s dumo. 
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Purpose: 
Deck Name: 
valling Sequence: 


Inout/Output: 


Subroutine PRINTR 


Prints a raw data file. 
PRANTR 


CALL PRINTR(STN,NSTN ,NTYPZ, ANNUM ,XJDATA , NP) 


[ FORTRAN 
1/0 Name Dimension Description 


STN 


NSTN 


NTYPs 


NNUM 


XJDATA 


Subroutines Required: 


Functions Required: 


Approximate Deck Length: 


* Dummy dimension 


Raw data array to be printed 


”. dicates station from which 
data was received 


Indicates data type 


Total number of words per raw data 
file 

Julian date (zero-hour U.T.) cor~ 
responding to the first data point 
within she file, 


Number of points per raw data file. 


None 
None 


260, = 178 
10 
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Subroutine CHRZNT 


Purpose; Sequences data chronologically in 
seconds from the reference date 
XDREF. 
Deck Name: CHRT 
Calling Sequence: CALL CHR@NT (STN,MAXSTN,NP ,XDREF ,XJDATA) 
Input /Output : 


Dimensions Description 


Data array containing a 
single file of raw data. 
STN(1,J) designates the time 
associated with the J-th point. 


I Dimension of STN. 
IL Number of points within the 

raw data file being processed. 
I | XDREF $f — weenee Reference date used for sequen- 


cing time, Value is 50 days 
less than XJDATA of the first 
raw data file processed, (This 
relatively short time interval 
assures that no loss of signifi- 
cance will occur when times 
differ by more than one second 
of time) 


XJDATA [| wenn Julian date (zero hour U.T.) of 
the first point within the raw 
data file being processed, 
Subroutine Required: None 

Functions Required: None 


Approximate Deck Uangth: 123, = 83,6 
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Restrictions: Time interval between successive raw data points 
must be less than 460 sec, 


Error Indicators: If time in seconds exceeds seven significant fig- 
ures, a diagnostic is printed and execution 
continues. 

Method: 


Prior to entering SUBR@UTINE CHR@NT, the chronological order of the 
data is mod 86400 seconds. The smoothing routine requires that each raw data 
file be chronologically monotonic (the mod sequencing would introduce 
discontinuities), and the sorting routine requires that the time for al] raw 
data files be with respect to the same date. These requirements arn met by 
selecting a reference Julian date (XDREF) 50 days prior to the date of the 
first raw data file processed and adjusting the chronological ordering of all 
raw data files to seconds from XDR&F. 


It should be noted that although XDREF is completely arbitrary with 
respect to the analysis (if a succeeding raw data file should happen to contain 
data with times prior to XDREF, the resulting negative time values would not 
affect either the smoothing or sorting), the computer system constrains the 
number of significant figures. Consequently, if time exceeds seven significant 


figures (9,999,999 sec), a diagnostic message, stating that XDREF may have to be 
increased, is printed. 
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2 SUBROUTINE ‘CHRONT (STN, MAXS TN, MB XDREF, XTDA TA) 


ENTER 


/o 


7G? s7 (4) 


IS 
OTREF *(XOREF~ XTDATA) ¥86400.0 


~|comeure TIME FROM REFERENCE DATE 
_|To ZERO HOUR THIS FILE. 


R22 
ie FOR EACH RAW OATA 


rT BOSE Eee POINT WUTHIAL THIS FILE. 


Sayer 25 4 a COMPUTE TIME INTERVAL BETWEEN 
loz ASS[STA (4 2) RESEMT AND PREVIOUS POINT: 


Ce a oe IF DT >86000,A ZERO HOUR 
LIES WITHIN) THE INTERVAL OT 
DELT=0ELT +86400.0||THUS TIME MUST BE INCREASED 
| BY ONE DAY (DELT). THIS TEST 
35 ASSUMES THAT THE ACTUAL TIME 


STA(1,1) = STN (1, 1) soeer| INTERVAL BETWEEA POINTS IS LESS 
_ITHAN HOO SEC. 


7B = STN(1, 2) 


sr (.t) +07 R«F | COMPUTE AND STORE TIME 
= STN(1, 1) tOTREF IN SECONDS FROM XDREF 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
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{ Subroutine FIT 

Purpose: Sets up the raw data for the smoothing routines. Also 
includes an option for discarding or adjusting raw data 
: not satisfying a specified criteria. FIT is predicated 


on the assumption that the data are recorded at intervals 

of approximately one second. If the time interval is much, 
much larger than one second, FIT should be reviewed and 
fewer than the presently considered 20 data points processed. 


4329 
Si aerhaxwneedlinne: 


Deck Name: FITT 


soma OS 


Calling Sequence:CALL FIT (STN,NP,IFI,AA,NPAA ,MAXSTN ,MAXAA ,NTYPE,NOAVG) 
Input/Output : 


Dimension Description 


ce ee EN 


er 


6 ,MAXSTN aba array containing a single 
file of raw data. STN(J,I) 
refers to the J-th element of 
the I-th data point. 


J code. 


- ener 


> 


moe 


1, time in seconds 
2, X coordinate * 
3, Y coordinate 

4, Z coordinate 

5, not used 

5, not used 


I NP Number of data points within 
the STN array. 


* X, Y, Z refer to as many as three general pieces of observed data. A3 
used in the differential corrections program, these data can be: 


1) Range 
2) Doppler 
3) Azimuth, elevation 
4 4) Range, doppler 
5) Range, azimuth, elevation 
6) Doppler, azimuth, elevation 
r 
as, ~6h,- 


SID 65-1203- 2 


~ = ean ET A 

Pe Win pM eee ne he 
Note A FS a RON a eae eae, 
ae a 


- [) Re |S ee ee Boel ot ” 


FORTRAN 
Dimension 


Deseription 


Number of coortinates per point 
within the SIN array. 


Data array containing the 
smoothed voints corresponding 
to a single raw data file. 
AA(J,I) refers to the J-th 
element of the I-th smooth 
data point. J code same as 
STN above. 


NPAA Number of smoothed points in AA 


MAXSTN Dimension of STN 


ener Dimension of AA 


—_——- Indicates type of data, See 
SUBRGUTINE PRGCES and footnote. 


Error flag. Non-zero value 
indicates that an initial valid 
point was not determined. 


Subroutines Required: SM@QTH (smooths a raw data segment) 
Functions Required: None 
Approximate Deck Lengths 1374, = 76446 


Nomenclature! 


FORTRAN 
Name 


Description 


Transfers a segment of valid raw 
data points from the STN array to 
the smoothing routines. torage 
1s allotted bv: 
D(NT), time coordinate - first 
valid raw data point 
D(NT+1), time coordinate ~ 
second point 


° 
e « 


D(NX), X coordinate - first 
point. 
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Dimension 


ERRSIG 


TADJST 


II 


IZVER 


MY MZ 


NT ,NX,IY ,NZ 


NOMPRE 


XAVG 


$ -~66— 


Description 


D(NX+1), X coordinate ~ second 
point. 


D(NY), Y coordinate - first 
point 


D(NZ), Z coordinate ~ first 
point 


Tolerance used with adjustment 
option. 


Adjustment option flag. Non-zero 
value mechanizes adjustment 
ootion. 


D array counter 


krror indicator mechanized if 

the number of smoothed data point 
exceeds the dimension of the AA 
array (NAXAA). 


AA array counter 


Index, within STN, of point being 
tested. 


Index, within STN, of last valid 
point tested. 


index, within STN, of coordinate 
to be adjusted. 


Index, within STN, of remaining 
coordinates, 


Index, within D array, of initial 
storage location of T,X,Y,2. 


Flac, Non-zero value indicates 
no more raw data in file. 


Averaze value of 40 coordinates 
immediately following coordinate 
being tested. Used in determini 
the first valid point. 


SID 65~1203- 2 


ask 


FORTR! HN 
Name Dimension Description 


Value of valid coordinate, 


Value of coordinate being 
tested. 


Error Indicators: if, using the adjustment option, an initia) valid point 


c2n not be determined using the averaging procedure, 
control is transferred back to SUBRAUTINE PROCES, an 
error message is printed, and the next raw data file 
is read into the program, 


If the number of smoothed points exceeds the dimensions 
of the AA array (NPAA> MAXAA), the first MAXAA points 
are loaded, an appropriate diagnostic is printed, and 
execution continues. 


Method: 

Prior to calling SUBRGUTINE FIT, a file of raw data observations 
from 4 single tracking station has been read into the STN array and has been 
chronologically adjusted (seconds from the reference date XDREF). 
This routine then partitions the raw data into segments and, after optional 
preliminary aata adjustments, smooths e«ch segment. When all raw data within 
the file has been smoothed, control is returned to SUBRGUTINF PRYCZS with the 
smoothed data in the AA array, 


ixamination of the raw data utilized in the sample problem disclosed 
random irregularities of azimuth observations heing recorded 369 degrees out 
of phase, especially prevalent for angles approaching 360 degrees (e.g. ,+354° 
recorded as ~ 6°), Consequently, an option has been included which adjusts these 
values before smoothing. A single raw observation is compared with the last 
valid point determined, and, if the result is not within a specified tolerance 
(ERRSIG), the observation being tested is adjusted by 360 degrees. The adjusted 
value is again compared and, if still not within tolerance, is considered invalid 
and discarded, The next raw observation is then tested and the procedure is 
repeated until the above criteria have been met, The point is then considered 
valid and stored in the D array. After the required number of points have been 
loaded into D, the raw data segment is smoothed and the resulting point is stored 
in AA. This procedure is repeated until the ruw data file has been exhausted, 
The initial valid point required to start the adjusting operation ig determined 
by sequentially comparing an observation with the average value of the 40 points 
immediately succeeding it until the tolerance criteria is met. 


Controis within the routine are set for azimuth adjustment. By 
resetting NX equal to the appropriete index within the STN array, and adjusting 
the tolerance Csrnstc) if’ required,any of the observed coordinates may be 
adjusted. To demechanize the adjustment option, set IADJST=0, 
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Since the time interval between raw data observaticns is assumed to 
be of the order of one second, the 20 points smoothed corresponds to a 
very short time interval. Consequently, if data are not collected in this 
mode, smoothing criteria should be adjusted. 


~68~ 


SID 65~1203- 2 


NYINLIS 


TILSEN AY FF O12 oo00s 


FALVMIAAOOD FHL AO 
‘NLS MWIHLIM ‘NX FQN 
SHL St XW Kvaaey 
GJ FHL NI SLINIOd 
GWA 10 FOVYOLS 
any inIWssneay 

HLNWIZY wos 
STOWLNOD 13S 


OFF? 


O59 


(TC) NLS = ('chvy 


ae 
O27 
O#LNWNQ7 


eee | 
ft IAT = IT4f 
INIHLOOWS 
LOOHLIM KILIIAIG WL ol 


MYY VOT ‘few F/ VL¥O NLOOWS LOW O@ 


Lfs INILNOUANS 


arrives} red 


FIdTI ‘t=L C49 LO Je! 
ee 
( advV7 =r O8990 — 


FNNILNID PP — 


peer 


Loman os 


-69-- 


SID 65~-1203- 2 


| Ee Pee cee 
Z 

4 

“4 SUBROUTINE FIT (Cour) 


THE FIRST YALIO RAW 
; DATA POT 18 DETERMINED 
WITHIN) THIS SECTION 


10 
: KGUNT» KOUNT +1 “] SET UP CONTROLS FOR 
K: Kel COMPUTING THE AVERACE 
& XSUM=0.0 VALUE OF THE ¥O POINTS 
3 < KEND =K+39 To Bh Tasrard. 
| 5 If AVERAGING REQUIRES 
4 KR 2 ’ MORE POINTS THAN ARE AVAILABLE, 
ra) KENOD ANP xes | KEND = NP READJUST NUMBER OF POINTS 
2 Be ee 70 GE AVERAGE D. 
< 50 : 
e DIVSR = KEND -Kel 
3 an ae IF NO POINTS ARE 
S YES AVAILABLE FOR 
= AVERAGING, RETURK 
é = ANO READ NEXT FILE 
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Subroutine SM@@TH 


Purpose: Smooths data to obtain the best fit (in the sense of 
unweighted least squares) to the polynomial 
Y=A+BX+CX@ 

Deck Name: SNGg 

Calling Sequence: CALL SM@OTH(X,Y,N, CYNST) 

Input/Cutput: 


Dimensions Description 


20 Observed data 


20 Values of the independent variable 
(time) corresponding to the observed 
Ys 


Number of elements within X and Y¥ 


Unknown constants to be determined. 
CENST(1) > A 
CANST(2) —~ B 
cgNsT(3) —~C 


Subroutines Required: MATMPY (matrix multiplication) 
MTINV (matrix inverse) 
TRANSP (matrix transposition) 


Functions Required: None 
Approximate Deck Length: 4739 = 315.5 
Restrictions: ND 3 
Formulation: 

The problem to be addressed is as follows: 

"A series of J observations (J>I1,I=degrees of freedom of function) has been 


made-over some interval of time (or another independent variable) of a linear 
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function in the unknowns of the problem, and an optimum estimate of the un- 
known parameters in the least squares sense is to be constructed,” From this 
statement of the problem 


Ytme +7 (2) 
where 
A 
¥ = e = J-vector of observations 
ti 
4 = J by I array of numbers relating the unknown parameters of 


the problem and the observations. 


x 
Xa 
X ede = I-vector of unknown parameters 
3 
Xx; 
—_ —_— ate 
0 : Y TRUK ~ \cesates 
Z- Mx 


At this point in the development, a scalar function of 7 must be constructed 
which can serve as the comparison function for various "curve fits" and which, 
when minimized will yield the "optimum" solution in the desired sense. This 
function, in its simplest form, is 

J 

Tr 
F2 Y) Ay = 9" (2) 

is ¢i 7 2 

However, since the various observations of y may be made with different 


precision, it is desirable to provide the capability for weighting the data 
inversely as the square of the variance in the observation, 


i.e, . 
ol 
pe (SE. 2V 
where 
Tye* © © J 
-l. ya 
V oOo “ +O | = symmetric 
(o) co) % | 


This second form will be utilized in the discussions which follow, 
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Substituting from the definition of 7 into that for F yields 


F{y - Mi} VV" fg - Mi 


> ut ote’ eal af ote T < oh 
z gv YX mary - Pym ex M VOMX 


Which when partially differentiated with respect to ®* (or ¥”) and equated to 
zero yields 


on- Y vv Me Xk MM 

— =i x al -! 
Xs? Y V'MEM VM) Ms 
ise " a 3 
Xwas =(m™v"M) "mly ly 


This solution, of course, contains the case where the data are unweighted 
(i.e.,v" = the identity matrix). 


Yi. 2(m™) MTY (4) 


This formulation, applied to the raw data, performs a vrelininary smoothing 

which assures more rapid convergence of the differential corrections solution 

in the main program and simultaneously reduces the amount of data to be processed. 
This step is accomplished by noting that the true tiajectory is nearly elliptical 
(i.e. a second order curve) and may be represented to a good degree, over a short 
time interval, by a parabola of unknown coefficients. i.e., 


ura +bx tox" 


! ¥, ¥ 
! x x a 

er fi ) (BPE OL} 
b o¥,0 8 


Consequently, if a series of observations (say ever a 20 second segment of the 
trajectory ) are utilized as successive values of y; and if the corresponding 
values of bi are utilized as the comnonents of xX , the unweighted least squares 
solution for the constants {a} may be obtained by equation (4). However, it 

is noted that since the values oft; may be large and differ in one of the final 
digits (e.g., data recorded in universal time at one second intervals could 

have values of & like 86300,86301,°"*) the columns of M will appear to be nearly 
linear dependent. This fact could result in severe numerical problems. Thus, 
it is suggested that the values of Xi be equal to t;- Emmdpomé (where Emigpomé 

is the value of 4, at the tenth point in the segment begin processed) to assure 
the maximum degree of significance in the computations. 
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Subroutine MATMPY 
(matrix multiplication) 


Purpose: MATMPY is designed to multiply any two conformable 
single precision matrices (with less than 70 elements 
each) to obtain the single precision product. All 
operations interior to the routine are performed 
in double precision to control roundoff and loss of 
significance. (The size limitation for this routine 
is a direct result of the establishment of several 
temporary storage arrays in the routine. The logic 
will admit any size matrices). 


Deck Name: MXPY 


Calling Sequence: CALL MATMPY (C, I, K, D, K, J, CD) 


Input/Output: 
FORTRAN | Math Common/ 
Name Name Dimension | Argument Definition 
I (rows) Arg Array of numbers 
K (columns) to be used in the 
premultiplication 
of D by C 
I D D K (rows) Arg jArray of numbers 
J (columns) to be premulti- 
plied by the matrix 
C 
I cD cD I (rows) Arg {Product array 


J (columns) | 
Subroutines Required: None 
Functions Required: None 


Approximate Deck 
Length: 1307, = T1lio 
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Purpose: 


Deck Name: 
Calling Sequence: 


Input/Output: 


mn net et 


Subroutine MTINV 
(matrix inverse) 


MTINV computes the inverse of a nonsingular square array 
(of up to 36 elements) using the theorem which states 
that if a sequency of row operations will reduce a 
matrix to the identity matrix, then the same series of 
cperations performed on the identity matrix will pro- 
duce the inverse. (Error control is maintained internal 
to the routine with double precision arithmetic though 
input and output are single precision. This double 
precision capability is directly responsible for the 
restrictions on the size of the matrices. No logic 
restrictions are involved). 


INV 


CALL MTINV (B, ES, N) 


Definition 


The N X N array of 
numbers to be inverted. 


(single precision) 


The N X N inverse of 
B (single precision) 


| The dimension of B 


as 


Subroutines Required: CHCOSE (Check for singular B) 


Functions Required: 


Approximate Deck 
Length: 


None 


The = 48470 
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Subroutine CHOOSE 


Purpose: CHOOSE is tilized in conjunction with MTINV to 
determine if the matrix (of order less than 6) is 
sufficiently non-singular to allow the inverse to be 
constructed without excessive numerical difficulty. 


Deck Name: CHSE 


Calling Sequence: CALL CHOOSE (A, E, M, N) 


Input/Output: 
aan Math “Common/ 
Name Dimension Argument Definition 


NHXN The array of sunbers 
which is being 
reduced to the iden- 
tity matrix 


E NXN Arg : The inverse being 
| constructed from A 


- ak Arg A row counter for the 
operations being 
performed on A 


- 1 Arg Order of the square 
array A 


Subroutines Required: None 


Functions Required; None 
Approxiias ave Deck 
Length: 271, = 185, 5 
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% Subroutine TRANSP 
(matrix transposition) 


array of numbers 
Deck Name: TRSP 
Calling Sequence: CALL TRANSP (A, N, M, B) 


Input/Output: 


i N (rows) 
M (columns ) 


M (rows) 
N (columns ) 


CS 


Subroutines Required: None 


Functions Required: None 
Approximate Deck 
Length: 137, = 9559 
‘ ~106~ 


Purpose; TRANSP constructs the transpose of an arbitrary 


Definition 


Array to be 
transposed 


The array containing 
the transpose 
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Subroutine GSPRT 


Purpose: Sorts the smoothed data. Primary ordering is time, 
secondary ordering is by station. 


Deck Name: NAOI 
Calling Sequence: CALL GS@RT (NA,N,M,K,L) 


Input/Output: 


1/0 Name Tescription 


Yo NA First element of data array 
to be sorted. 


T N Number of groups (points) 


I M Number of elements in each 
group (point) 


I K lement number in the group on 
which to sort. K*M 


I L Number of passes. Positive 
value, sort from maxinun to 
Minimum value. Negative value, 
sort from minimum to maximum 
value. For a complete sort of 
N items, L=N-1. 


Subroutines Requir i: None 
Functions Required: None 


Approximate Deck Length: 367, = 247, 
a 


8 0 
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Restrictions: 


The data is considered to be a one dimensional array 
composed of multiple, equal lengthed data blocks, each 
block consisting of consecutive elements. Consequently, 
if a specified element is to be sorted in a multi- 
dimensioned array, the proper ordering of the indices 
is unique. 

e.g. Assume the second element is to be ordered ina 

2 by N array, then the input array must be of the form 
A(2,N}, not A(N,2) 
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Subroutine CHANGE 


Purpese: Transforms units of the smoothed and ordered output 
data 

Deck Name: CHANG 

Calling Sequence: CALL CHANGE (B,MSTART ,MADSUM,DELT ,C,g) 

Input /Output : 


FORTRAN 
Name 


Dimensions Description 
Data array to be converted. 
Cooresponds to either STN 

or A array in SUBRSUTINE PRECES 


Location of first point in B 
array to be converted. 


Losation of final point in B 
array to b3 converted, 


Difference between program 
reference date and reference 
date used for sequencing time. 


Coafficients required for coord, 
transformations (Km, Ki/sec, rad) 


Output data array containing 
smeothed, sorted, and transformed 
date. Cooresponds to a single 
logical record on the output 
tape to be read by the main 
program. (FS4-507) 


a Pt A 


er 


eee ae 


* Dummy dimensioned 
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Subroutines Required: 


Functions Required; 


Approximate Deck Length: 


Error Messages: 


ANGM@D 


CXA 


CXAPLB 


656. = 430 
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None 


2195- 


(transforms time from sec to 
days) 


(insures that angular measurements 
are positive, mod 360 deg) 


(transformation of form D=O-A. 
sed for converting azimuth and 
elevation from degrees to radians) 


(transformation of form D=-C« A+B. 
Used for converting doppler 
reading to range rate.) 
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= 7O TW, TF. 
30 4 
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DECK LENGTH IN GCTAL [S$ 
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Subroutine TWTF 


Purpose: Converts chronological time (sec) as referenced to 
XDREF into days relative to 1950.0 (JD 2433282.423) 
and divides the resulting time into whole and frac~ 
tional parts, 


Deck Name: TWTFR 
Calling Sequence: CALL TWIF (7. DELT, TW, TF) 


Input/Output: 


FGRTRAN 
Name 


Time in seconds from reference date 
Used in chronological data sequencing 
(XDREF in SUBR@UTINE PRGCES) 


Dimensions 


Diffe.ence betwean program reference 
date anc reference date used for 
sequencing data, i.e, DELT=XDREF~ 
XJDREF 


Integer number of days from 1950.0 
(JD 2433282,423) 


Fractional number of days defining 
epoch 


Subroutines Required: BESSEL (Besselian calendar correction) 


Functions Required: None 
Approximate Deck Length: 14, = 7815 


Error Messages: None 
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Subroutine BESSEL 


Purpose: Special purpose routine which converts time from a 
reference date of JD 2433281,5 to that relative to 
JD 2433282.423 

Deck Name: BES 

Calling Sequence; CALL BESSEL (TW, TF) 


Input Output: 


FORTRAN 
1/0 Name Dimensions Description 


Input: integer XJDREF, (XJDREF = 
Julian Date Jan. 0, 1950, 
2433281,5) 


Output: integer days from 1950.0, 
(1950.0 corresponds to Julian Date 
2433282423) 


Fractional day 


‘ Subroutines Required: None 
Functions Required: None 
ba 
Approximate Deck Length: 55, - 45.9 
Error Mess2ges? None 
uaa 
~126- 
SID 65-1203- 9 
} 
ania seal hon ROOTS 2 Lotter Repro 
ra “ sheen, 
fe RE 


00¢0C0S39 
06T00S38 
08 t00s3g 
OL100S38 
09100S38 
ostucssag 
0+ 100838 
O¢ 100S39 
octooss?2 
OTTOOS38 
00TO0GS38 
96000838 
08000S38 
02000838 
09000S38 
05000S38 
07000839 
0£000S39 
occoos3g 


S8/82/60 


~ 


ONS 

NuNnLay 

€26°O - si = dl 

ot O41 BS 

220°O + di = Si 

“tT - AL = AL 

S20! 99 (€¢6°0°39°H1L)d1 
(4i‘AL)73SS38 SNILNOVANS 


E2v°7BZEEV? “AF OL SGNDdSAIVUND O°OSET NVIWSS39 
S°teceeve = O°0S6T AiVG NYIINEG = sayvgrx 
*3°3 
*YVGNS TVD NVIVSSS3E SHL HLTA GNOdS3uyOD OL =i GNV AL 
SLD3UYND BNIANOY SIHL *O°OS6T 31LVO 39N38333U WVYUINUd 3BHL YO 


*3SDduNd 


**«e T3SS30 JNJANDYENS «ee 


~- (S3NJ1I - $34 


ysocsd 


AN3BW31VLiS 308NDS N330 ~ 


xaxe VSOC-¢SA +42 


N 


VLIGOUOLOUWO © 


went 
ve 


ae: A 


ee 


ww Ser " 
Pe EMT ARR IIRR - E cot en te mse ey 
= ae 


rw 
er 3 ee 


hoe 


> 
Mee 
ar 


“oe 


a 


: 


rs 


, 


Does 


‘ % ok 


. 


a 


3 


2 


N 
t 
: 
t 
IN 
Oo 
i=) 
Hi 
wn 
46000 SI W198 NI HLONSI W330 
22000 ve O€ €Z000 v9 oz 02000 VS Si 
$1000 VY OI 42000 Vi SZ 01000 v1 S 
NOILVIO1 Nil N33 NOILVIO4 NJI N33 NO1LVIOF Nal N33 
3INIGNOdS3 YAOI Nd] N33 
ae : 4 NOIL935 DOISAS 
G371V2 S3NELNOYENS 
€ NO1193S 13S$S38 
WT et Ree, OY oe We” 6 = SAiNIDd AYING 
; 73SS3@ 3NILnoOyNaNs 
© 8 
sks dVW 39VUOLS $39 = m 
~ a 68/82/60 VSOESS | 
4 


PE ee 
5 abe om: & 2 


? ate gee ; 


3 i te a . Eo SA ts b : nes . oe ats 
a ee. 


Function ANGMZD 


Purpose: Tnsures that angular measurements are 
positive, mod 360 deg.. i.e. 0°<B<360° 


Deck Name: AMD 
Calling Sequence: B=ANGM@D(A) 
Input/Output : 


Description 
Angle to be made mod 360 


Anod360 a 


Subroutines Required: None 


Functions Required: None 

Approximate Deck Length: 75, = 61.5 

Error Messages: None 

Restrictiona: ~ 3600 €A< 3960 (Range of A may be increased by 


increasing the upper index limit of the Dg loop) 


#% Amod360 = A < by where ( ) represents the integral part of the 
360 


quotient. 
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Function CXA 


Purpose; Transformation of form D=Cra. Used for converting 
azimuth and elevation from any arbitrary measure. 

Deck Name: CA 

Calling Sequence: D=CXA(C,A) 

Input/Output: 


$$. 


Description 


Variable to be transformed 
Conversion coefficient 


Cra 


Subroutines Required: None 


Functions Required: None 
‘ Approximate Deck Length: 33, = 2745 
Error Messages: None 
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Function CXAPLB 


Purpose: Transformation of form D=C*AtB. Used for converting 
doppler reading to range rate (Km/sec). 
Deck Name; CAB 
Calling Sequence: D=CXAPLB(C,A,B) 
é Input/Output : 
FORTRAN 
1/0 Name Dimension Description 
I A Variable to be converted 
x B Conversion coefficient 
Z CG Conversion coefficient 
H Q D D=C- A+B 
ii 
A 
4 Subroutines Required: None 
a 
a Functions Required: None 
a Approximate Deck Length: 36, = 3045 
‘ Error Messages: None 
i 
4 
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SUBROUTINE PRINTS 


Purpose: Prints the logical record of smoothed and ordered 
data being written on the output magnetic tape. 
Deck Name; PRNTS 
Calling Sequence: CALL PRINTS (%,NPERGP ,XJDREF) 
Input/Output : 
FORTRAN | 
Name Dimensions Descriptions 


1/o 


71% Array containing the smoothed 
and sorted data to be printed. 
G(J,I) refers to the J-th 
| component of the I-th data 
point. 
J index code: 
1, integer days from J.D, 
2433282,423 
2, fractional day (U.T.) 
3, station name indicator 
4, data type indicator 
5; X coordinate *% 
6, YX coordinate 
7, 4 coordinate 


NPERGP Number of points within the 
@ array. 
I/0 XJDREF Program reference Julian date. 


Subroutine Required: None 


Functions Required: None 


% Dummy Dimension: See SUBRGUTINE PRGCES 


wt X,Y, Z are dummy variables used to describe the observed data (not 
rectangular coordinates) 
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Restrictions: 


Approximate Deck Length: 


Advantage was taken of the fact that the logical 
record size basing written on the output magnetic 
tape happens to correspond to a printed page. 
Consequently, if XM@DS is changed in SUBR@GUTINE 
PROCES, headings will not necessarily be printed 
in the proper order, 


224, = 148), 
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SUBROUTINE PRINTS (G,NPERGP, XIDREF) 
ENTER | 


DATA arene a 7 
GERD INITIALIZE LINE COUNTER 
| Caan eg cepa eat iS ace owed 
Ss vf 
HEADING 
NEW PAGE, PRINT HEADING 
/04 
po 0d 35 I+ 1,NPERGP a DO FOR EACH POINT 
j 
| IS 
| _] WWCREMENT LINE COUNT 
| 
| 20, 
| . 
BD (3,7):NSTN FIX FLOATING VARIABLES 
| B(4,Z):nTyYP 
30), 
LINE, TW, TE. NSTH, N x, 
. Po tead DR ReIGL a | PRINT A DATA POINT 
ete J 
| 28 
L.————{ COW TINVE 
40 
R e 
Resa ETURN AFTER PRINTING 
A SINGLE LOGICAL RECORD 
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PROGRAM OPERATION 


This section describes the general requirements for program operation, 
For a more detailed discussion, see the sample problem (Appendix B). 


Input Data 


Input data to the program are entered on FORTRAN 10 digit decimal data 
sheets, Although other data sheets are acceptable, the latter is suggested 
because its format precludes the possibilities otf input errors. (See Figure 2) 


Columns 2-3 and 4-5 of the first data card are reserved for inter indica- 
tors which control the printing of raw (PR) and smoothes (PS) data respec- 
tively. To demechanize the print options, insert zero, The first colum of 
this card must contain a blank unless this is also the final card (i.e. using 
preset data exclusively), 


The remaining cards, excluding the final card, are utilized for inputting 
the following information (one station per card). 


Description 


Blank (indisates data follows on this card) 


Integer identifying the station for which the 
remaining data on the card applies, Same value 
as NSTN on the input data tape. 14 K # 10 


STNAME(K) Order in which the Kth station is to be sorted. 
This number must be consistent with the orderin 


employed in the differential corrections program. 


Range conversion (R=C,R; to Km) 


Range rate conversion (R=CoR + C34 to Km/sec ) 
Azimuth and elevation conversion (A=C,A; to rad) 


The final card contains an integer (non-zero, non-blank) in column 1. The 
remaining fields are blank, This card indicates that all data have been 


entered, 
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It should ba noted that, although the first and final card are necessary 
for normal operation, only those data cards for which stations are specified ao 
are required, Further, the ebove locations are all presst within the program 
(see SUBRGUTINE PRESET), I2% the preset value ia desired, the field is left 
blank, 


Input Taps 


The program requires that tracking station observations be input on a mag- 
netic tape with the following format and specifications (see Figure 3), ' 


l. FORTRAN IV, binary mode, multi-file tape. 
2. Data arranged one physicai file per station per pass. 
3, Two logicai records per file, 


A. First logical record - 4 word information record containing the 


following: 
NSTN indicator which identifies the station from which data 
was received, 1< NSTN € 10 
NTYPE indicates data type * 
i=: range 


2= range rate 

3= azimuth, elevation 

4= range, range rate 

Se: range, azimuth, elevation 

6-: range rate, azimuth, elevation 


NNUM total number of words per file (including information 
record} 


XJDATA Julian Date (zero hour U.T.) corresponding to the first 
piece of data within second logical record, The lead-~ 
ing characters "24" have been omitted from all Julian 
Dates for the sake of numerical significance, 


B, Second logicai record - logically packed tims and coordinate data. 
Each raw data file must be chronologically sequenced, However, 
the sequence dozs not have to be monotonic increasing, ¢.g., the 
observations from a single station for a single pass (single raw 
data file) may be ordered in universal time 84901, 84902, ... 
86400, ly 2, 3 aes 


4, Input tape is mounted on logical tape drive 8, 


5, End of tape is indicated by an information record containing zero's 
immediately after the final physical file containing data. 
% Jf agimuth and elevation are recorded, the coordinate system must be 
topocentric with respect to the plumb line. 
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RECGRO 
ONE STATION 
PER FILE PER DATA 
PASS RECORD 


Figure 3 


INPUT TAPE FORMAT 
RAW DATA. 


NTYPE 


TRACKING STATION 

TYPE OF DATA 

NUMBER OF WORDS PER FILE 
JULIAN DATE OF FIRST DATA POINT 
TIME 


COGRDINATES DATA 
AS INDICATED POINT 
BY NTYPE 


NSTN7O SAIDICATES NO MORE 
BATA ON THIS TAFE. 


° FORTRAN & 

e BINARY MODE 

* MULTI- PYSSICAL FILES 
"LOGICALLY PACKEO Tiwk 
ANOQ COORDINATE DATA 


SID 65~1203- 2 


Variable Dimensions 


Although auxiliary, intermediate tape storage is available, maximum pro- 
gram efficiency and optimum core storage utilization may be achieved by proper 
adjustment of certain key array dimensions in the MAIN routine. To facilitate 
the determination of the optimum set of dimensions, the principle storage 
functions of these arrays are summarized below. 


Auxiliary Tape 
Storage Not 
Mechanized Remarks 


Auxiliary Tape 
Storage 
Mechanized 


Array Name 


raw data 
(single file) 


raw data 
(single file) 
smoothed data 
(all files) 


smoothed and 
ordered data (all 
files) 


STN(6 ,MAKS IN) 


A(6,MAXA) smoothed data smoothed data loaded from AA 
| (at least single (all files) (file by file) 
file) 


smoothed and ordered | if the number of 

data (all files) smoothed points 
exceeds MAXA, 
auxiliary tape 
storage is mech- 
anized by the pro- 
gram 


if auxiliary 

tape is mechanized, 
transfers data 
from tape to STN 
array for ordering 


AA( 4 MAXAA) smoothed data smoothed data transfers data 
(single file) (single file} from smoothing 


routines to A. 
(file by file) 


If an array should be underdimensioned, during execution the program will 
usually stop loading when the array is filled, compute the proper dimension, 
print an appropriate diagnostic, and either continue execution without loading 
the remaining data or terminate execution with a core dump. 


= Pe 
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Storage Limitations 


Assuming that optimum core storage has been allocated by proper adjust- 
ment of the STN, A, and AA arrays, there are two upper constraints on the 


1, 


260 


amount of data thay may be processed. 


The amount of raw data per station per pass (one physical file on the 
input data tape) exceeds the capacity of the STN array, This problem 
is easily corrected by partitioning the data into smaller files, i.e. 
the input data tape is arranged one physical file per fractional pass. 


The amount of smoothed data exceeds the capacity of the STN array. 
This is highly unlikely but could occur if an excessive number of 
physical files are read in at one time, In the program's present 
form, smoothed data are read into the STN array from temporary tape 
storage (the auxiliary tape mode would have been automatically mech- 
anized) until the number of smoothed points equals the dimension of 
the STN array, The required dimension for STN is then computed (it 
may be possible to readjust dimensions) ard a diagnostic is printed, 
The data in STN are then ordered by the sorting routines and execu- 
tion continues as normal, The excess data still on the tape are not 
utilized. 


Tha latter problem may be corrected by: 


i 


ae 


Da 


decreasing the amount of raw data per computer run. Two disadvan- 
tages of this method of operation are: it may not be advisable to 
separate the raw data if adjacent stations have overlapping spheres 
of influence, or it may be necessary to arrange the raw data files in 
chronological order on the input tape. For example, since each run's 
data would be sorted separately, the chronological ordering of the 
final output tape may net be correct, 


decreasing the number of smoothed points by increasing the ratio of 
raw to smoothed points, Presently, a segment of 20 raw data points 
is smoothed to produce a single smoothed point. 


It is anticipated that future versions of this program will include an option 
where a filtering process would selectively discard the surplus number of 
smoothed data points, 
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PROGRAM OUTPUT 


Primary output consists of a magnetic tape which is read by SUBROUTINE 
DATAPE in the program (FS4-305). Optional printout of the raw and/or 
smoothed data is also available (non-zero values for PR and PS) as explained 
in the PROGRAM OPERATION section. 


Output tape 


Mtilined below is a summary of the format and specifications of the 
output tape (See Figure 4). 


1. FORTRAN IV, binary mode, single file. 


2. Data arranged by logical records. 
A. First logical record - 4 word information record containing the 


following. 
XJDREF 


NGPS 


NPERGP 


NPREM 


UD 33282.423 (Jan 0 1950) 


Total number of logical records on the tape (excludes 
information record). 


Number of points*® per logical record (excluding final 
data record). 


Number of points in final record. 


B. Remaining records - Data records containing, 


NTY?PE 


ep Fu7 


Time. Integer days from 1950.0 or JD 2433282,422 
Time. Fractional day. (U. T.) 

Station from which data was received. 

Type of data. See ITYPE in SUBROUTINE DATAPE. 


Components of the observation vector. 


* A point is defined to be the ordered set of words TW,TF,NSIN ,NTYFS,X,Y,zZ 
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FIGURE 4 


FS4-305A 
OUTPUT TAPE FORMAT 
SMOOTHED AND ORDERED DATA 


XJDREF 


PROGRAM REFERENCE DATE 


INFORMATION | NGPS | TOTAL NUMBER OF RECORDS 
"| NUMBER OF POINTS PER RICOKD 
RECORO NPERGP |“ UMGtuniNa FivaL RECORD 


NUMBER OF POINTS IN SINAL RECORD 
WHOLE PAYS FROM 1950.0 
FRACTIONAL) JO. 24393292. VAI 
STATION 


TYPE OF DATA DATA POINT 


AS INDICATED 


DATA RECORD CONTAINING: GY NTPs 


¢NPERGP POINTS 
* SEVEN WORDS PER POINT 


COORDINATES ; 


° FORTRAN WW 
* BINARY MODE 


APPENDIX A 


The enclosed routine (SUBROUTINE DATAPE) reads the magnetic tape 
generated by the preprocessor progran FS4~305A. Details concerning the 
magnetic tapes format and specifications may be found in the PROGRAM OUTPUT 
section of the FS4-305A write-up. This routine is an integral portion of 
the operation of the differential corrections program FS4-305 and is also 
documented along with the routines of the Data Filter Group of that 
document. 
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Subroutine DATAPE 


Purpose: Reads a specially generated magnetic tape containing 
smoothed and ordered coordinate data, 

Deck Name: DAPE 

Calling Sequence: CALL DATAPE (TW,TF,ISTN, TYPE, DATA, KGUNT) 

Input/Output: 


Description 


Time. Integer days from 1950.0 
(JD 2433282,423) 


Time. Fractional days, 


Station from which data was 
received, 


Indicates type of data in @DATA. 
Code: * 

1, Kange 

2, Range Rate 

3, Azimuth, Hlevation 

4, Range, Range rate 

5, Range, Azimuth, Elevation 

6, Range rate, Azimuth, 
Elevation 


Coordinate data. 


Control indicator. 
Code: 
1, Do not return next point 
(INPUT) 
2, Return next point (INPUT) 
3, No more data pcints on 
tape (OUTPUT) 


%* If azimuth and elevation are recorded, the coordinate system must be 
topocentric with respect to the plumb line. 
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Subroutines Required: None 
Functions Required: None 
Approximate Storage Required: 664g 


Restrictions: Requires that the input data tape generated by 
FS4-305A be mounted on logical tape 9, 


Nomenclature: 
FORTRAN 
Name Dimensions Description 


7,36 Buffer array containing a single 
logical tape record. 
First vass indicator 


Logical record counter, 
IGP=1 ,NGPS 


Number of logical data records on 
tape. 

Number of points per logical record, 
excluding final record, 

Number of points per final record. 


Buffer array point index, NPAINT=1, | 
(NPERGP or NPRM) 


Program reference Julian date, #* 


Method: 

The first CALL statement to this subroutine mechanizes the input data tape. 
An information record containing the progrem reference date, number of logical 
records, and number of data points* per record is read into the program. The 
first logical data record is then read into the wffer array A and, after 
extracting the first data point from the buffer, control is returned to the 
calling routine. Subsequent CALL statements extract single points sequentially 
from the buffer, After the final point within the buffer has been extracted, 
the next legical record is read in and the procedure is repeated until all 
data has been read from the tape. 


*A data point is defined as the ordered set TW,TF,ISTN,ITYPE /PDATA. 
HID (24) 33282.423 


-155-~ 
SID 65-1203- 2 


Soa Reda cl ord “a ee poy Te ne pam ee ERE 


wee 


ee apna RET: dota Nan ine f 


Areay wrsttT 
onns THOIFY vsvd tSuleS avs 
(6) 207+ 


sé 


"GHOITH 
rents Ml SLN/Od 
yO MTITHOAN Of Vv007 ¥ 
oyorsy HIS FLNIOd 


é FgvL 
wo PHOITY VLVI 
INO KINL THOW 


FO YIQWAN LTS | of 
aquoosy (6) 
aly Ya#d SLNIOd yo YPFIWAL © gz OMIM Fe Sf oN 


‘gx0n Te THNIS QNITATIXT 
‘oyoosy UPA SPYOM IO ysewnn 
sovorszey vLvo 70 yszawan Wt ole 
FLVG BSINIAILTIA wry D0td « 
CHOI NOL eOAN[ COLT wnos ovir 
(6) FAVL 


IAVLVG INILNOYGNS 


@ 


YFLNNOCD LN100e ansesae onl | 


os 


s7A 
wav, VLVO 


—— 


onw eanscoows WOUyr avsy 
778 WoITeY LSWS SVH 


of 


2 ENO vLVg 
LXIN NAALIY 


s 


wIhWHNOD 7 
wILNng? 


ztyos 


WoDF ov ® 
LNIId @ 
Lsvlse 


-156- 


SID 65-1203- 2 


“PHOIST oY 

WAS WI SLAIOCS FO 
ywoewanv os WwAnZ 
THOWIY WIS FLIMOS 
SO WIFPSWON LIT 


SLNICH FJYOW ON | 
—~ YOLVICH LIS 


$9 


WILNNOD LNIOd LISTHY 
06 


Avouwv WISING 
OLW/ THOITH vivagv dvIH 
(&) FoVL 
cy) 


ON 
ae éEgvoD7ey TWN/S 


SIA WL FG StHL WM 


08 


wILNWNOD 
THOIse LNINFYIN/ | 


SL 


ZS 


é 7¢7hL VO 


ON WSLNIOd 280W ANY ~~ O09 


ON 


(LNOD) FdVLYO INILNOYGNS 


NYOAL TY ) 


SE 


SILVNIAMWOOD 
VLE SO FAL @ 
NOWILWLE o 
IW1IL 8 

Avaewl wIAS/AY wows 
LMI TIONS V LIVYLXT 


. GO 


STA 


ékvyHV wz77I/NF W/ 
~ SLNIGS FUCW ANY 


GF 


@ 


-157~ 


SID 65-1203-2 


wh ree ae erecta enter a bie ne nc 


ee 


ee I AN N — 


a Se ee ee Sein meee a7 


_ __ 02€00000 eA) ; ee serene ye aD ee 
09€00000 “6 LINN SAIUG Jdvi WoOISTT NO G3IINNOW Ce) 
0S€06000 39 VSOE-%S3 AG AZIVYINID Jdvt Viva INdNI SHI LVHL SI¥INDIY ae 
rea 3400000 we 
Q€€00000 SSNOTA9DTYES3IY 3 
—“OZ€00000 7° = — ———“— Ss ae Page ee eee ag een Oe eg CRE ake QE Pah 
01T€00000 “3dVIVG TW) HOVA AB Y3ASNG SHL WOU is) 
——-———_Huea0000 ~~ OFLBVYIXI-ATININGSS NSW JIV SINYHd VIVO FTINIS “Y AVYIY OD 
06200000 Y333NG@ BHI OLNY 3dVi WOVS GV3Y Juv SGHOIIY Viva WIIS08T =. 
892700000, i Ap aT Poe 
oz2z0c0000 *QO0H13W 3 
———————~“G9200000,.—~—~—~—~—‘“TCSC*:*S pi eee ey eee eg a Oa Ge et ee ae Sere ee sy 
05200000 CINdINOI°AdVI NO SINIOd VLVG SYOW ON °E 9 
942000007 andy” CIN “*XAN NUAIY 7C ae 8 
0€Z00000 CINdNIT} *INIOd LXIN NUfLay LON OC ST 4 
Oe2e00000 ~~ =" "SFA TODIGNY TOYUINDGD FINNS Y ~~” 4 
01200060 *¥viVG JIVNIGHDOHD ‘vivao 9 
5200000, NOLLVAD 132 *HiNWizv*’sivd JONVT 9 = 7 7 4 
06TC00GO NOTLVASTS*HINWIZVESONVY °S 9 
——_——— “45100060 7 8 SSN eSSNVE TS ve : 4 
OL1TO00GOO _  _LNOTAVAZTSSHIANTZY SE 9 
99100000  °»~— ; ——" “3T9u Fonvy #200 4 
06100000 cath : 3SONVY *T 9 
——sey ———“GeTo0000 OTT eLVOO NI ViVa 5O 3dAi SZRVIIGNT S3dALT ss) 
0€ 100000 *G3ASIIIY SVM VIVO HOTHM WHYS NOTIVIS * NISI 3 
OZ TOOOOCO AVG WNOTLIVYUS “AWIE & St ) 
OTTO3dvO *3LVG JONIVIIZY WYYNOYd NOUS SAVA YAIZINI °3HIL 6 ML 9 
Te eT TOSdWOT ho ee tig ee ee ioe cia — 9 
06003 dvd *SNOFZINIS3Q INJWNOUV it; 
“= " O8003d 9d : : Fete So ere a ye ean 6] 
0L003d¥vG “viv J3IYNiOYOOS G3IyITYO GNV 3 
- -"  * QHQQOZdvVa CSHLOONS SNINTVINOD JdVT JTLSNOVW GZivYINID ATIWId3adS V SAV3y 9 
0so03dva ten aticte: eterna ty ethan 0 7 ; 2 
——"—"""“959g0adVG — | Fasoddnid i) 
o<cooadvG Ee re as . i) 
~  “9z003adva0 eae JADVIVG ANTIINOYERS #4 wee SOE—>S4 #249 ‘ | 
eae 2 A Sales eed ‘s Oh sete vice _. es . +3 
wo 
ee (USINAT = INBW3LVES a9unes N37 - __. 36¥0 Rat Break. ' 
S3/9T/E6 SO€+Sa 
4 
#3 
ae 
nai Sh ett SS ee oe , ‘ t E 


pee > | IC Ne, Sa ne a 


+a 


04200000 cone ee oe HANTOUNE ZIV = _SL OC1 
0€ 206060 CANIDAINS ET) = Mi S6 
___0¢L00000 = °¥BasdING HONS ANTOd 29VMIXS ne eee + 
01200060 I = £NIOGN 06. 
T¥ 00200000 . Si CONSNST=IS G24 tare cr ervey) (63 aYaY SB 
06900000 W3UdN = GNIN € SJONSDR*dO! DAT 08 
__. 98900000 hi8.73 ; ; ; ; 14 d9l = d9?I G2 
62900000 <= 
_ 09900000 - eee ct oe GET 91 09 OL 
06900000 eres yee ey Tg, =< € = ENNOY $9 
____. 0%900000 See sete dts ek ~ at ; SZ O81 09 € SAHIN*SN°dOT FSF 09 
C€900000 — © ~~~" “" *NuNsay GNV 


02900600 € 982i tWND3 INNA L3S *Viva 2u0W 

~OT900000 ON 3¥ “*dYHIIY WwII90T IXSN v3" cane! 

00900000 G6 O01 99 ( GNIN*’STPINIDIN JAY SS 
~ VESOONGN ~O9F9 IXSN AV3Y *YIHING “SHL WOHT a 


O8SOOOOODSLIVULXZ NIIG SVH INTOdD TWNIS at 


““GLS00000 ~~ " 1 + ANIOdN = INIOHAN OS 
09500000 *INNOD INTOd YISING ENIKIYONI 
0$500000 1 = d5I Sy 
LY 0%S00000 _ (ONIN THI ee tart ers yve} (6) GV3Y OF 
~OESOONCD eee : d9SYIGN = GNSN SE 
0zS00000 WIUdN = dIYIIN (T°DI°SdDINIAT OF 
CT OT SOOODG TT RITIN STSITIN SAINT S3SUGLX «6(6Y AVI SZ 
00500000 1 = 2SYIST OZ 
I W6y00000 = Se nae metal mira rine One? eevee a SE GUN Mad Gt 
08+00000 *OupI3Y VAVO ISYI4 ONV 
——"GT 550000 “NOTIVREDINY DVIS "SSVI LSU 
09%00000 Of 92 09 ¢ OP IN°ESUT4T IAT OT 
— OS900000 - SP oT pd T TsseiNnOyY Fars” 
0400000 
PEDO ee ena. Se SOS ee ea 
02400000 / Ow€ / dOY *INTOAN “ISUTST VIVG 
=O 1SOOO0G ry eg ee ut, tye ee ag 
0000000 {€)vivao * (9€42)v NOTSNSWIG 
Pr OOGEOU00N ra ee ee, Sem a ee ae : 
____oB€00000 = LANMOM EV EVODSIGALIANISTASL‘MLIZdVIVG BNTINOWNS 
— (S3NST — SNBW3iVIS 39uNOS NSS -  S3dvO 
“68/91/80 ©” © GOE#S3 


YY YO YVuUY 


Bare tee Te 


N 
5 
' 
uN 
wa 
a 
Lape] 
se 
028000006 ie 7 ad ONE Viet ogee!) fe ee ee ne 
01800000 : NUNL3Y OFT 
~ O0806000 ~ ~~ IEE EEE TTT SNPINDD SET 
0620000008 CANTODNSZIV = CEPVAYGO S21 
08200000 ~~ CINTOGN 98 = (e)vivan Oet 
02200000 ANTON ASHY = CTPVIVOD STI 
09200000 : = ~ CINTOEN* SIV = 3dkLT OTT 
os200000 NT ON ED = NIST SO% 
oS 
Moe SINST O~- ENGWBEVES BDUNOS NZD - ada 
S8/9T/80 SO0ESS4 


t 
t 
S1soo0 W62 Ss C1S00 vlz So 19400 ss ii OF 
bs$00 vot S€ 25400 vel o€ SE400 v6 D7 sz 
tO or et, ee ees OOO NE Sie, een 8 Uie BNGOO 8 nt NEE oe OS 
£2400 vs OT 12900 v99 Set £1400 vi a S 
NOTAVIOT = =—-_ssNGT _ .N#3 NOTLIVIOD NST _. _Nd3 NOFEVIOT OO ONS NAS, 
ee. . JONZGNOGS3YINI NST NAB Siar cone oe Pena oa cena 
~%t ~~ NO1ISSS ~~ ———OCSD ISAS €t  NOVIDSS5 337007 
et UNOTEOBS EPID NOTEDIS OPI OT NOT NDS 
6 NGT1995 *1093° 2 NOT15395 “183° 2 NOTi93$ “wu 
_9 So NO@THIFS CC GONAE OCF GO NOTEDIS BS a. NOTEISS OF AMUST 
O371V2 S3NYENOUENS : ee 
ee e”~C«NOLDSS)C”~<“«é‘i! 
TT GENYOd AVN i sts—<“—sSs—C(‘; fae we 
I B07 00! a ANIOON yo cs &, SOOO a ES cee ete et ORO te 
a ey 0. T61)) ONAN I 10%00 W3UdN T 60400 d5USAN 
els eee eS ae COD a ON es ete ie eee Oe eee os ee ee ge aE 
3441 NOTIVIO7 WEWAS 3dAL NOTLVIO7 TOEWAS JdAk NOTIVIO TOSWAS 
i a eas a SSTGVIUVA WVUIOUd G3NOISN3WIGNN .. ——r—CC—C— oy oye 
Sa Ba RS seg aS he Eg ihe Gk a Pp a i aca a ee COU ited rea at 
IdAi NOTLW901 TWEWAS FdAl NOT1VI01 TOEWAS adi NOTLVID? TDOWAS 
S3TAVIUVA HVUNONd GINOISHSWIG See ee ee ia a ea 
= a is BONANO  ANESOD SONS oo, Steet ate | led ee le als 
ree gO ON SONS Fe it yt eh es oe 
$8/91780 SOEeSS 


OE Tn a OT 8 
aaa ts ke 

bast 3 
tach 


“Bo” 


s1p|65~12034 2 
& 


ae Bee, Bese eee ee wa 79900 SE W390 NE HLONDT 4930 
vi9 OvT LT900 v49 SZ1 St¥900 _veg _02zT 
¥09 Sit 90900 ves OTT 10960 V9 ——~ 661 ; 
V5S oot 02500 VIS 06 14S00 vty ss 8B? ‘ 
veo 08 0€S00 VoE oy] 92500 Wot ——" $9 i 
VLE GL __ 22500 vzE 09 Z72S500 __ ves _ «SB it 
N 7 
dvH FOVUOIS __ 3dv0—“(‘( ie i, 
$8/91780 SOESS : 
Pee! 
ae aces kT Cs ‘ : eo [ et sd * : 
é 


APPENDIX B 


SAMPLE PROBLEM 


Input Tape 


Input data for the sample problem consists of two naper tapes 
(generated by a Packard-Bell 250 Computer) containing the following observa- 
tions for the ECHO II satellite passes 6069 and 6070: elevation and 
azimuth in degrees, doppler reading in kilo-cycles, a lock-on indicator, and 
universal time in seconds. Since the IBM 7094 computing system utilized 
does nct have a paper tape input capability, a short IBM 1401 program was 
written to transfer the data directly to magnetic tape with the format 
shown in Figure3 (one physical file per pass). The leading information 
record for each data file contains the following indicators. 


PASS PASS 
SYMBOL DESCRIPTION 6069 6070 
NS IN Code indicating the tracking 1 2 
station recording the data 
NTYPE Code indicating the type of data 6 6 
(range rete, azimuth, elevation) 
NNUM Total number of words per file. 3564 3392 
This was determined by the 1401 
program. 
XJDATA Julian date (zero hour U.T.) of 38850.5 38850.5 


the first observation for each 
file. (arbitrarily selected for 
the sample pro’-lem. ) 


A third file containing NSTN equal to zero indicates to the program that 
811 data has been read from the tape. 


Graphical representation of the raw data is presented in Figures 1-B 
through 6-B. The same data has been printed out in the sample problem 
output (raw data files 1] and 2). Examination of Figure 5-B will disclose 
randon. irregularities of azimuth observations being recorded 360 degrees out 
of phase, especially prevalent for angles approaching 360 degrees (e.g., +354° 
recorded as -6°.) These points are adjusted in subroutine FIT immediately 
prior to smoothing the data segment. 
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This figure also serves to illustrate a data format which is not 
acceptable to the program. Although the data is generally sequenced 
chronologically, a few random errors in time values were generated by the 
tracking stations (e.g., raw data point No. 245, raw data file No. 2, 
sample program output, Page ys 


It is the user's responsibility to insure that the time data are recorded 
and sequenced properly. The tracking errors in the second raw data file of the 
sample problem were not corrected. The smoothing routine (see Subroutine FIT) 
processes segments of 20 consecutive points (time is not explicitly utilized) 
and assigns the time value associated with the midpoint of the segments to 
the resulting smoothed point. Since the erroneous time values did not occur 
at segment midpoints, the program was not affected. 


Execution time for the sample problem was 1 minute and 18 seconds. 
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Input Data 


Since the preset values were utilized exclusively in the sample 
problem, ne input data was required. This sample serves to demonstrate the 
operation of the preprocessor. However, the data tape generated is not 
suiteble for direct application to the differential corrections program 
(FS4~305) because a unit conversion was employed in conjunction with 
the doppler data, and because of the assumed date of the passes. The 
deck setup consisted of the $ DATA control card followed immediately with 
a "6" card (non-zero value in column 1 indicates end of data.) However, 
for illustrative purposes, Figure 7-B has been included to show how the 
data deck could have been set up to achieve identical results. 


Variable Dimensions 


The required values for the variable dimensions were determined by the 
procedure outlined t-low. (See "Variable Dimensions" and "Storage 
Limitations" in the "Program Operation" section.) 


AA(4,MAXAA): This array must be large enough to hold the smoothed points 
corresponding to a single raw data file. The largest raw data file contains 
3564 words. The number of words per information record is 4, and there 

are 4 words per data point (time, doppler reading, azimuth, elevation.) 
Thus, the number of raw data points per file is (3564-4)/4 = 890. Now 

the smoothing routine reduces 20 raw data points to a single smoothed point. 
Therefore, since 890/20845, MAXAA was dimensioned by 50. 


A(6,MAXA): The primary purpose of this array is to receive the smoothed 

data from the AA array. Since MAXAA was set to 50 and there are two files 

of data, MAXA was set to 2 X 50 = 100. If there had been storage limitations, 
MAXA could have been reduced (tne smoothed points would then be temporarily 
stored on tape.) However, this array must be large enough to hold the 
smoothed points corresponding to a single raw dats file. 


STN (6,MAXSTIN): The usual criterion for dimensioning this array is 
number of raw data points per file. The sample problem has MAXSTN = 1000. 
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Output 


Primary output is the magnetic tape to be input into the main program 
(FS4-305). Printed output includes the following. 


Input Data; 
The input data is printed out on the first page. Note that the 
printout for the enclosed sample problem contains the values preset by 
the program. ' 


Raw Data: 


By setting the raw data print indicator equal to a non-zero value, 
each raw data file will be printed out immediately prior to being smoothed. 


Smoothed Data: 
By setting the smoothed data print indicator equal to a non-zero 


value, each logical record of the output data will be printed immediately 
after being written on the output tape. 
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